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In recent years, a close connection between the description of open quantum systems, the input-output for¬ 
malism of quantum optics, and continuous matrix product states in quantum field theory has been established. 

So far, however, this connection has not been extended to the condensed-matter context. In this work, we 
substantially develop further and apply a machinery of continuous matrix product states (cMPS) to perform 
tomography of transport experiments. We first present an extension of the tomographic possibilities of cMPS 
by showing that reconstruction schemes do not need to be based on low-order correlation functions only, but 
also on low-order counting probabilities. We show that fermionic quantum transport settings can be formulated 
within the cMPS framework. This allows us to present a reconstruction scheme based on the measurement of 
low-order correlation functions that provides access to quantities that are not directly measurable with present 
technology. Emblematic examples are high-order correlations functions and waiting times distributions (WTD). 

The latter are of particular interest since they offer insights into short-time scale physics. We demonstrate the 
functioning of the method with actual data, opening up the way to accessing WTD within the quantum regime. 


I. INTRODUCTION 

Continuous matrix product states (cMPS) have recently 
been recognised as powerful and versatile descriptions of 
certain one-dimensional quantum field states [1-3]. As 
continuum limits of the matrix product states (MPS)— 
a well-established type of tensor network states underlying the 
density-matrix renormalisation group machinery [4] —they in¬ 
troduce the intuition developed in quantum lattice models to 
the realm of quantum fields, offering similar conceptual and 
numerical tools. In the cMPS framework, interacting quantum 
fields such as those described by Lieb-Liniger models have 
been studied, both in theory [1, 5, 6] and in the context of 
experiments with ultra-cold atoms [7]. 

On a formal level, continuous matrix product states are in¬ 
tricately related to Markovian open quantum systems [1, 2]: 
The open quantum system takes the role of an ancillary sys¬ 
tem in a sequential preparation picture of cMPS. Elaborating 
on this formal analogy, cMPS can capture properties of fields 
that are coupled to a finite dimensional open quantum sys¬ 
tem. This connection has been fleshed out already in the de¬ 
scription of light emitted from cavities in cavity-QED [2, 8] in 
the quantum optical context, under the keyword of the input- 
output formalism [9]. 

Another methodological ingredient to this work is that 
cMPS have been identified as tools to perform efficient quan¬ 
tum state tomography of quantum field systems [7, 10-13], re¬ 
lated to other approaches of tensor network quantum tomog¬ 
raphy [11, 14]. These efforts are in line with the emerging 
mindset that for quantum many-body and quantum field states, 
tomography and state reconstruction only make sense within 
a certain statistical model or a variational class of states. Im¬ 
portantly, in our context at hand, it turns out that cMPS can 
be reconstructed from the knowledge of low-order correlation 
functions alone. This is a very attractive feature of cMPS: 
Recently, a reconstruction scheme has successfully been ap¬ 
plied to data on quantum fields obtained with ultra-cold Bose 
gases [7]. Read in the mindset of open quantum systems, 
cMPS tomography can be interpreted as open system tomog¬ 


raphy by monitoring the environment of the open quantum 
system (see Section II). 

In this work, these methodological components will be put 
into a different physical context and substantially developed 
further as illustrated in Eig. 1. At the heart of the analysis 
is a tomographic approach, applied to an open quantum sys¬ 
tem, yet brought to a new level. In Section III, we extend 
the set of tomographic methods within the cMPS framework, 
showing that the dynamics of the ancillary system and of the 
whole open quantum system is not only accessible from low- 
order correlation functions, but also from low-order counting 
statistics. Specifically, we prove that for generic systems, the 
two density functions Pq and Pi—which express the prob¬ 
ability of detecting zero and one particle, respectively, as a 
function of the time since the last detection—provide suffi¬ 
cient knowledge to successfully perform tomography of the 
open quantum system.The physical application of the estab¬ 
lished methods will also be different from the cavity-QED or 
the quantum field context: Here, we treat fermionic quantum 
transport experiments within the cMPS framework. 

In a general transport setting, a scatterer is coupled to a 
left reservoir (the “source”) and a right reservoir (the “drain”). 
Eermions (with or without a spin degree of freedom) can be 
seen as jumping in and out of the scattering region from the 
source to the drain and can be described by a leaking-out 
fermionic quantum field. In Section IV, we show how the 
dynamics of the open quantum system (scatterer and leaking- 
out fermionic field) can be encoded into a cMPS state vector. 
To provide the reader with an intuition about the equivalence 
between the cMPS language and a more traditional Hamilto¬ 
nian formulation, we will consider one of the simplest setups 
in quantum transport: a single-level quantum dot weakly cou¬ 
pled to two reservoirs. These results are also valid for trans¬ 
port experiments of ultra-cold fermions between a “hot” and 
a “cold” reservoir as recently realised in Refs. [15, 16]. 

This will clear the way for making use of the tomographic 
possibilities offered by the cMPS formalism to access various 
quantities in quantum transport that are not yet measurable 
with current experimental technologies (see Eig. 1). Emblem- 
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FIG. 1. Extension and applications of cMPS based tomography (see Sections III-V). Previous works [7, 10] have shown that measurements 
of low-order correlation functions Cn with n = 2,3 are sufficient to access higher-order correlations using cMPS tomography based on the 
cMPS matrices M and D (right side, black arrows). The first achievement of this work is the formal demonstration that measurements of low- 
order counting probabilities constitute an alternative to measurements of low-order correlation functions for carrying out cMPS tomography 
(blue arrow). We show that the probabilities to detect zero or one particle (Po and Pi respectively) are sufficient to reconstruct alternative 
cMPS parameter matrices M and P, from which higher-order correlation functions can be computed. The second achievement of this work 
is to extend the applicability of this cMPS machinery to quantum transport experiments. As a main illustration, we show that cMPS-based 
tomography provides an access to the distribution of waiting times W. The according statistics are not directly measurable due to experimental 
limitations on single-particle detectors. However, we demonstrate with experimental data that they can in fact be reconstructed from the 
knowledge of low-order correlation functions (broad blue arrows). 


atic examples are higher-order charge correlation functions 
and the distribution of waiting times (WTD, see Section V). 

The waiting time is defined as the time interval between 
the arrivals of two consecutive electrons. Therefore, the WTD 
provides a privileged access to short-time physics, short-range 
interactions and the statistics of the particles. As such, it has 
gained a lot of attention recently [17-24], but WTDs suffer 
from their difficulty to be measured effectively: Measuring 
WTDs requires the detection of single events while ensuring 
that no events have been missed—for instance, due to the dead 
time of the detector. 

With present technologies, WTDs in transport experiments 
can be measured when the injection rate of electrons is within 
the kHz range. Here, the current trace is resolved in time and 
the WTD can be directly deduced from it. This is for instance 
the case in the experiments [25, 26]. It is important to keep in 
mind that experiments in the kHZ range can only mimic trans¬ 
port properties of classical particles—no quantum effects such 
as the statistics of the injected particles or the coherence prop¬ 
erties of the scatterer can be observed at these frequencies. 
To observe quantum mechanical effects within the setup, one 
needs to move to the GHz regime, which can be achieved ei¬ 
ther with DC sources with a typical bias of tens of meV, or 
with periodically driven sources at GHz frequencies [27-31]. 
In the GHz range, the current trace can not be resolved in time 
so that the measurement of the WTD is not feasible at present. 
In contrast, second- and third-order correlation functions have 
been proven to be feasible, but are accompanied by exceed¬ 
ingly difficult measurement prescriptions [31, 32]. 


In the light of this discussion, we propose in Section V an 
indirect way to access the WTD with methods that are within 
reach of the experimental state of the art. Namely, the dynam¬ 
ics of the full open quantum system is accessed from measure¬ 
ments of low-order correlation functions (typically second- or 
third-order). This is made possible with a cMPS formulation 
of the transport experiments as explained in the following sec¬ 
tion. 


We illustrate this indirect path of accessing the WTD by 
considering real data obtained in the experiment of Ref. [26] 
where single electrons tunnel through a single-level quantum 
dot in the kHz regime. Both the current trace resolved in 
time and the two- and three-point correlation functions have 
been measured. Although no quantum effects are present as 
explained above, the data allows us to demonstrate a very 
good agreement between the WTD deduced directly from 
the current trace and the WTD obtained via our reconstruc¬ 
tion scheme based on the data of the correlation functions. 
This gives substance to our protocol based on cMPS to ac¬ 
cess the WTD with present technologies. We claim that this 
method remains valid in the GHz frequency range and for 
more complex systems such as a double quantum-dot cou¬ 
pled to two reservoirs—which would exhibit quantum co¬ 
herence effects—and for quantum transport experiments with 
fermionic quantum gases. 
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II. CMPS TOMOGRAPHY 

In order to present a self-contained analysis, we start by 
reviewing the cMPS formulation of capturing a finite dimen¬ 
sional open quantum system [2] and the tomography proce¬ 
dure of reconstructing the relevant cMPS parameter matri¬ 
ces [10]. Consider an open quantum system (in cMPS terms 
the ancillary system) with dimension d (called bond dimen¬ 
sion in that context) and interacting with one or more quan¬ 
tum fields that are described by field operators for differ¬ 
ent fields a. Its dynamics can in general be represented by 
different mathematical objects: 

a) The master equation in Lindblad form, which governs 
the evolution of the ancillary system described by its 
state vector |^) defined on the Hilbert space H of di¬ 
mension d X d. The degrees of freedom of the coupled 
fields are traced out in this approach. 

b) The set of n-point correlation functions of the coupled 
fields. 

c) The full counting statistics of the field system, i.e. the 
complete set of cumulants of the probability distribu¬ 
tion of transferred particles. The n-th cumulant of the 
generating function is linked to the n moments of this 
distribution, which correspond to the n-point correla¬ 
tion function. 

d) The cMPS state vector IV^cMPs), which we now intro¬ 
duce. 


A. cMPS reconstruction from correlation functions 

An intuitive way of establishing the cMPS state vector 
IV^cMPs) consists in starting from the well-known Lindblad 
equation. This equation describes the evolution of the state 
p in time via the Liouvillian superoperator £, 


p = ^[p] 

= -l[K,p]-\i^ {{RiRa,p] - 2Re.pRi) . 

a=l 

The first term relates to the free evolution via a Hamiltonian 
K G while the last two terms describe the coupling to 

the environment (the according operator is known as the dis- 
sipator). The matrices G a = 1,... ,p, correspond 

to jump operators between the system and external quantum 
fields The matrices K and {Ra} completely charac¬ 

terise the evolution of the system. 

Making use of the Choi-Jamiolkowski isomorphism [33] 
(which maps linear superoperators from Ri io R 2 to linear 
operators acting on Hi ( 8 ) H 2 ) the state p is mapped to a state 
vector Ip) and the Liouvillian C to the matrix T [1, 2] with 

T = Q*0l + l0Q + ^i?* (2) 


The matrix T G is known as the transfer matrix and 

the matrix Q is defined as 

Q = (3) 

a 

Formally, the isomorphism introduced above is defined by the 
following relations for an operator and the product of opera¬ 
tors 


p \p) 

, ( 4 ) 

pB ^ {A* ® B)\p). 

Being closely connected to K and {i?a} introduced above, the 
knowledge of the matrix and T and its components provides 
access to the dynamics of the open quantum system, and al¬ 
lows to directly derive the according Lindblad equation. 

The (translationally invariant) cMPS state vector |T^cMPs) 
on the interval [0, L] is defined in terms of the matrices 
Qi{Ra} and the field operators by 


iV^cMPs) = Tra: 


V exp J dx (g) i + ^ Ra 


( 8 ) 


1^). 

( 5 ) 


This expression is related to the path ordered exponential that 
arises when integrating the Lindblad equation. The embed¬ 
ding of the cMPS state vector IV^cmps) into Fock space be¬ 
comes clear when expanding the path ordered exponential 
V exp. For more details, we refer to Ref. [3] where the au¬ 
thors formulate the cMPS in different representations such as 
the Fock space and a path integral formulation. After inte¬ 
gration, the ancillary system is traced out via Trane and the 
resulting term is applied to the vacuum state vector \ft) where 
1^2) = 0 for each a. 

Compared to the Lindblad equation, the main difference is 
that the degrees of freedom of the ancillary system are traced 
out such that its dynamics is mapped into the dynamics of the 
coupled quantum fields The evaluation of expectation 

values of field operators leads to expressions that only contain 
quantities from the ancillary system, and information about 
the ancillary system can be inferred from according field op¬ 
erator measurements. For the sake of clarity, we restrict our¬ 
selves to the case where a single coupled quantum field, de¬ 
noted as 7^/3, is measured. 

The density-like correlation functions of the measured 
quantum field 'fp then read 

Cn(x) = (V^cMPsI ^(^ 1 ) • • • iV^cMPs) , (6) 

where x := (xi,..., Xn) and n := According to the 

calculus of expectation values in the cMPS setting [3], insert¬ 
ing Eq. (5) into Eq. (6) in the thermodynamic limit L ^ oc 
leads to the expression 


C'n(x) = 


( 7 ) 


lim Tr 

L^oo 


g-D(-£/ Xn) Xn — l)jy^ 
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With D we denote the transfer matrix T —introduced in 
Eq. (2) —in its diagonal basis, 

D = X-^TX, ( 8 ) 

where the columns of X represent the eigenvectors of T. 
Analogously, the matrix M denotes in the diago¬ 

nal basis of T, 


M = X-\Rl0Rp)X. (9) 

Let us mention that the knowledge of X is in principle not 
necessary to reconstruct the matrices Q and R and hence the 
according Lindblad equation [10]. 

Specifically, the second- and third-order correlation func¬ 
tions take the form 


(f 

C 2 {x) = Tr [e^“Me^“M] = ^ (10) 

i=i 


and 


(f 

C3{x,x')= y] 

3,k=l 


XjX^Xk{x'-x) ^ 


with {Xj} being the eigenvalues of T. Due to the translation 
invariance of the system, we can set = 0. The tomographic 
possibilities of the cMPS formalism can be understood from 
Eqs. (10-11): If the products {Mi/cM/e jM^ 1 } are known, we 
can—using gauge arguments [ 11 ] —require each of the matrix 
elements {Mij} to be equal to one, which enables us to ac¬ 
cess each Mkj by dividing the appropriate terms: 


Mi,iMijM^-1 “ 


( 12 ) 


Both numerator and denominator appear as coefficients in C 3 
and can be determined with spectral estimation procedures. 
This means that in principle we just need to analyse a three- 
point function in order to obtain the building elements M and 
D of arbitrary-order correlation functions. 

This reconstruction scheme demonstrates the central role 
of the matrices M and D to derive the different equivalent ob¬ 
jects that describe the dynamics of an open quantum system: 
the Lindblad equation, the set of n-point correlation functions, 
the full counting statistics of the number of transferred parti¬ 
cles and the cMPS state vector. These matrices M and D can 
therefore be considered as the central quantities on which our 
reconstruction machinery is based; this is illustrated in Pig. 1. 


third-order [25, 26, 32]. A priori, this might render the recon¬ 
struction of the matrices M and D infeasible, but the work in 
Ref. [11] proved that this limitation can be circumvented by 
making use of the structure of the cMPS state vector combined 
with the thermodynamic limit. 

Por a given finite region / and a fixed bond dimension 
d, all expectation values can be computed from all corre¬ 
lation functions Cn(x) taking values in the finite range /, 
X = (xi,... ,Xn) C This contrasts with the situation 

of having access to correlation functions Cn (x) for arbitrary 
values of X = (xi,..., but for low n. Here, arbitrary val¬ 
ues X imply the thermodynamic limit, i.e. the finite region / 
tends to infinity. Then indeed, low order correlation functions 
(typically Ci, C 2 , C 3 ) are sufficient to reconstruct an arbitrary 
expectation value of an observable supported on /. 


III. CMPS RECONSTRUCTION FROM LOW-ORDER 
COUNTING PROBABILITIES 


In this section, we extend the central role played by the ma¬ 
trices M and D for tomographic purposes by showing that 
they (more precisely: their equivalents M and V) are also ac¬ 
cessible from low-number detector-click statistics, i.e. the idle 
time probability density function Pq and the density function 
Pi, which correspond to the detection of zero and one particle, 
respectively, within a certain time interval r. 

It is well-known that correlators and counting statistics are 
closely related. When assuming perfect detectors, the proba¬ 
bility to observe n events in the time interval between t and 
t + r is given [34] by the expression 

Pn{t,t + T)= (13) 

2 ^_ ■y'^m—n rt+r rt-\-r 

T / ^ 7 TT / dti • • • / d-tjYi CfYi(t\^ ^ 2 ? • • • 5 tm) 5 

{m-n)\ Jt Jt 

where the correlation function Cm has been introduced in 
Eq. ( 6 ). Por a translationally invariant system, we can with¬ 
out loss of generality set t = 0. Purthermore, when changing 
the integration bounds and performing the limit L —^ (X), we 
obtain 


Pnir) ■■= P„(0,0 + t) = 

ptl ptn-l 

dt2--- dtnCn{T,ti,t2,---,tn), (14) 
Jo 




B. Use of the thermodynamic limit 

Intuitively, it is clear that the reconstruction of the matrices 
M and D should gain precision by increasing the number of 
correlation functions Cn on which the reconstruction scheme 
is based. The same statement is valid when increasing the size 
of the set of available counting probabilities P^. But in gen¬ 
eral, experiments will only provide us measurements of low- 
order correlation functions, typically those of the second- and 


with Cn{T,ti, ...,tn)- = 

(15) 

with the canonical unit vector ei, the diagonal matrix V of 
(5*(g)l + l(8)(5 with basis transformation matrix Y, A4j := 
Y~^(Rj 0 Rj)Y, and Z := Y~^X, where X diagonalizes T 
as defined in Eq. ( 8 ). With Eqs. (13-15), the low-order count¬ 
ing probabilities Po(t) and Pi (r) within a cMPS formulation 
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are given by similar expressions to Eqs. (10) and (11), namely 


(f 


Po{T) = e{Z 

i=i 


d? 

Pi(.t)= y~^ zi,jMj,kZk,i 

j,k=l 


X 


(1 ^j,k) 


Qk'kP _ Qk'jP 





(16) 


(17) 


with {jij} being the diagonal values of V and {zj^k} being the 
elements of the matrix Z (with inverse =: (%,/c)). See 
Appendix A for details. 

As a first step, we can extract from Pq and Pi the coeffi¬ 
cients {zi^jZj^i} and the eigenvalues {/ii}, which give rise to 
V. The matrix elements of M. can then in principle be deter¬ 
mined using gauge arguments and under the assumption that 
the additive components of P^ are linearly independent. From 
M and P, the cMPS matrices Q, P, and K describing the dy¬ 
namics of the open quantum system can be determined in a 
straightforward way (see Appendix A for details). 

Let us comment on the feasibility of this reconstruction 
scheme with present technology. In order to measure Pq and 
Pi, efficient single-particle detectors without dark-counting 
and tiny dead-time are necessary. Dark-counting leads to de¬ 
tector output pulses in the absence of any incident photons 
while the dead-time is the time interval after a detection event 
during which the detector can not detect another particle. Al¬ 
though significant experimental efforts have been made in or¬ 
der to improve single-photon [35] and single-electron detec¬ 
tors [36, 37], the state-of-the-art for single-particle detection is 
not yet sufficient to perform a reliable measurement of Pi. For 
the moment, these experimental constraints make the recon¬ 
struction scheme based on Pi only valid on a formal, mathe¬ 
matical level. In the light of the recent experimental progress 
towards the reliable detection of single particles, we believe 
that this idea will become relevant in the future. 


IV. APPLICATION TO FERMIONIC QUANTUM 
TRANSPORT EXPERIMENTS 

Very recent works have successfully formulated experimen¬ 
tal setups in cavity QED and ultra-cold Bose gases as well 
as the corresponding measurements in terms of cMPS [7, 8]. 
This allowed them to make predictions for higher-order corre¬ 
lation functions that are not accessible experimentally and to 
investigate the ground-state entanglement. 

Here, we tackle the problem of formulating quantum trans¬ 
port experiments and the corresponding measurements (av¬ 
erage charge current, charge noise) in cMPS terms. To this 
end, we demonstrate that the field that is leaking out and is 
measured in a quantum transport experiment belongs to the 
cMPS variational class. We then provide an example to il¬ 
lustrate the equivalence between an Hamiltonian and a cMPS 
formulation by considering one of the simplest transport ex¬ 
periment, namely single electrons tunnelling through a single- 
level quantum dot. We derive the first-order and second-order 


correlation functions in cMPS terms, and show that we re¬ 
cover the well-known expression of the average current and 
charge noise, when writing the cMPS state Eq. (5) in terms of 
the parameters of the quantum system. 


A. Quantum transport experiments 
in terms of cMPS 

We now turn to a description of the physical setting un¬ 
der consideration. We assume here transport experiments, 
where single electrons transit through a scatterer coupled to 
fermionic reservoirs. The reservoirs, considered at equilib¬ 
rium, are characterised by their chemical potential and their 
temperature via the Fermi distribution. The bias energy and 
the bias temperature between the different reservoirs will set 
the direction of the charge current. For the sake of simplic¬ 
ity, we restrict ourselves to two reservoirs, the source and 
the drain. This transport setting can be described by a tight- 
binding Hamiltonian, 


— ^sys + ^res + ^int, (18) 

where ^sys relates to the quantum system under investigation, 
which acts as scatterer. It is characterized by discrete energy 
levels Si with occupation number operators given hy d\ 

(di^a- and dl ^ denote the fermionic annihilation and creation 
operators for an electron on the energy level i and spin de¬ 
gree of freedom a = t, i). The Hamiltonian relates to 
the left and right reservoirs, and Hint describes the interaction 
between the quantum system and the reservoirs, 

E.= E E / dEEcijE)c^,,{E), (19) 

a=L,i?cr=t4'^^ 

Pmi= E E / d-E A,t,a(-E’)rft,<T 0C^,c.(-E’) +h.C.) ■ 

a=L,R i,a 

( 20 ) 

The creation and annihilation operators of the reservoirs, Ca 
and , satisfy the canonical anti-commutation relations and 
a = L^R denotes the left and right reservoirs respectively. 
The amplitude ta,i,cr sets the interaction between the quantum 
system and its environments. 

In order to model a DC source, the energy levels in the 
left and right reservoirs are assumed to be densely filled up 
to the energies Ep + eV and Ep, respectively. Here, Ep is 
the Fermi energy and V is the bias potential applied on the 
’’source” reservoir. At zero temperature, the bias energy eV 
enables uni-directional transport of electrons between the left 
and right reservoirs. It plays a similar role to the frequency 
bandwidth when, e.g., considering cavity QED setups, and 
fixes the energy domain over which electronic transport takes 
place. 

With this assumption about the direction of propagation of 
the electrons (from left to right), we will see that Eq. (18) is 
equivalent to a generalised version of the cMPS Hamiltonian 
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introduced in Refs. [1,2], 

^cMPs = Q (8) i + (^Rl (8) '0L + ^ 5 (21) 

where the matrices Q and {Ra} and the quantum fields {^a} 
have been introduced in Section II. The cMPS Hamiltonian 
for quantum transport experiment reflects the direction of the 
current: a fermionic excitation present on the left of the scat- 
terer is annihilated at the scatterer as described by the quan¬ 
tum field (an electron jumps into the scatterer). Similarly, 
a fermionic excitation present on the right of the scatterer is 
created at the scatterer as described by the quantum field 
(an electron jumps out of the scatterer). The case of a multi¬ 
terminal setup can be considered in a similar way. Showing 
that Eqs. (18) and (21) are equivalent implicates that there is 
a fermionic quantum field leaking out of the scatterer to be 
measured and that it belongs to the cMPS variational class. 
Such a description of the transport experiment corresponds to 
a fermionic version of the input-output formalism of cavity- 
QED setups. 

Using Eq. (20), the quantum field leaking out of the quan¬ 
tum system, ^R(t), can be written in terms of the creation 
operator in the right reservoir cr; the incoming quantum field 
can be written in a similar way in terms of the creation opera¬ 
tor in the left reservoir cl, 

[ dE a = L,R. (22) 

JE 

The Eermi sea for the electrons is taken into account in the 
following way: On the right side of the scatterer, the quantum 
field satisfies ^^(t) |E^f) = 0, where \Ef) denotes the state 
of the Eermi sea at energy Ep, whereas on the left side of 
the scatterer, '0L(t) + eU) = 0, where the state vector 

\Ey + eV) defines the state of a Eermi sea at energy E^ + eV. 

Assuming that the energy levels Si of the quantum system 
are well inside the bias energy window eU, we can rewrite the 
integration over the energy domain E as dE = dE. 

This assumption is the so-called large-bias limit, which is 
considered in order to derive the master equation correspond¬ 
ing to the tight-binding Hamiltonian. In quantum optics, it 
corresponds to a finite frequency bandwidth, which allows 
the use of the rotating wave approximation [8, 9]. In the 
following, we assume that the interaction amplitude is spin- 
and energy-independent within the interval [E’f, Ey -\- eV] : 
ta,i,a{E) = ta- Let US remark that the demonstration remains 
valid with an interaction amplitude that depends on spin and 
energy. 

In a rotating frame with respect to the energies of the reser¬ 
voirs and after a Jordan-Wigner transformation using the def¬ 
initions of the quantum fields t/^r^l given in Eq. (22), the 
Hamiltonian in Eq. (18) can be rewritten as 

= Hsys (8) i + E E + h.c.)(23) 

q;=L,R i,cr 

Eollowing quantum optics calculations—which remain valid 
in this case because Ht is a transport version of the spin- 
boson model—we finally arrive at an effective non-Hermitian 


Hamiltonian 

= j ffsys - ^ ^ I ® 1 (24) 

\ a=L,R i,(j J 

i,a 

with ta •= Expressed in the eigenbasis of H^ys, the 

operators s/T^di^a- and ^ take the form of matrices 

labelled i?R,i,cr and i?L,i,cr, respectively. The effective non- 
Hermitian Hamiltonian can then be rewritten in a compact 
form 

H,ff = Q0i (25) 

+ E ® ® ^R,a W) ' 

i,a 

When comparing this effective Hamiltonian with Eq. (21), the 
identification of the matrix Q and the matrices {Ra} is di¬ 
rect. Eor spin-less fermions, the matrices R verify . ^ = 0 
in order to satisfy the Pauli principle. Eq. (25) demonstrates 
that transport settings can be adequately formulated within the 
cMPS framework. This result is important as it clears the way 
for applying methods from cMPS tomography to fermionic 
quantum transport experiments. 


B. Single energy-level quantum dot 


To illustrate the input-output formalism and the cMPS for¬ 
mulation of quantum transport experiments, we consider one 
of the simplest setups, namely a single energy-level quantum 
dot, without spin-degree of freedom, weakly coupled to two 
fermionic reservoirs. Even though this experiment is charac¬ 
terised by Markovian dynamics, this example is of particular 
interest for this work as it has been widely investigated ex¬ 
perimentally. In Section V, we will use real data obtained in 
Ref. [26] for this setup to show that cMPS tomography allows 
us to access the electronic distribution of waiting times. 

This simple transport experiment is sketched in Pig. 2 and 
the corresponding Hamiltonian reads 

HT=sd^d^ ^ f dE (^^yI^d(S) cl{E) E HrQS‘ 

a=L,R'^ 

(26) 

Because of the weak coupling assumption between the dot 
and the reservoirs, the tight-binding Hamiltonian Ht reduces 
here to a tunnelling Hamiltonian. Assuming that we perform 
a measurement on the right of the scatterer, the first two corre¬ 
lation functions of the right quantum field ^R(t) read in terms 
of cMPS matrices 

(4^r) = lim Tr 0 Rr)] (27) 

L^oo 

and 


(^R(0)^R('r)^R(r)V’R(0)) 


= lim Tr 

L —>-CXD 


.T(L-t) 


(i?^®i?R)e'"(i?^®i?R) 


(28) 




7 



FIG. 2. (Color online). Scheme of a transport experiment through a single-level quantum dot. a) The single-level quantum dot with energy £ is 
tunnel-coupled to two biased reservoirs with coupling strengths Fl and Fr. Spin-less single-electron tunnelling events take place in the energy 
window eV above the Fermi sea with energy Ef . b) The same transport experiment from the open quantum system perspective for a cMPS 
formulation. The single-level dot is described by Hsys and coupled to fermionic quantum fields and The coupling matrices T^l/r 
depend on the parameters Fl/r, see Eqs. (24-25)). The transport direction fixed by the biased energy between the left and right reservoirs is 
ensured in the cMPS formulation by imposing + eV) = 0 and ^r|F^f) = 0. 


The matrices T^r/l correspond to the operators ^/T^d and 
expressed in the eigenbasis of the single-level quan¬ 
tum dot, {|0), |1)} (empty and occupied state), 


Rl 


0 0 
x/Tl 0 


, = 


0 

0 0 J ' 


(29) 


Inserting these expressions into Eq. (27), we recover the 
well-known expression for the steady-state current of a single- 
level QD coupled to biased reservoirs [38, 39], 

(4^r) = =: (J)ss. (30) 

Furthermore, we can derive the noise spectrum from Eq. (28) 
via the MacDonald formula [40-42] 

This example aims at bridging the gap between a more tradi¬ 
tional Hamiltonian and the cMPS formulation, which allows 
to write these well-known expressions in terms of the param¬ 
eter matrices Q, T, and {Ra}- 


r —expresses the probability of having detected zero particles 
in the interval [0, r]. In terms of Pq, the WTD has first been 
derived in the context of quantum transport experiments in 
Ref. [19], 


W(t) = (r) 


d^Poir) 

dr‘^ 


(32) 


Here, (r) denotes the mean waiting time. Inserting Po{t) in 
cMPS terms (Eq. (16)), we arrive at an expression for W in 
terms of the cMPS matrices D, V and Z defined in Eqs. (8) 
and (15), 


WM = -el - 2DZ-^V + Zei. 

c 

(33) 


The normalisation factor c > 0 ensures that >V(r)dr = 1. 
Equation (33) allows us to access the WTD from the measure¬ 
ments of the low-order correlation functions via the use of the 
cMPS machinery to reconstruct the cMPS matrices D, V and 

Z. 


V. RECONSTRUCTION OF WAITING TIME STATISTICS 

In this section, we address the problem of accessing the dis¬ 
tribution of waiting times in electronic transport experiments. 
As mentioned in the introduction, a direct measurement of the 
WTD is not yet possible due to the lack of single-particle de¬ 
tectors with sufficient accuracy. Here, we propose to recon¬ 
struct the WTD based on the experimental measurements of 
low-order correlation functions. The reconstruction is carried 
out using the machinery based on cMPS presented in Sec¬ 
tion II and the formulation of transport experiments in terms 
of cMPS as exposed in Section IV. 

A. Definitions 

The statistics of waiting times can be expressed in terms of 
the probability density function Pq, which—as a function of 


B. Results based on experimental data 

We demonstrate our novel approach to derive the WTD 
from the measurement of correlation functions using exper¬ 
imental data obtained in Ref. [26] for spinless electrons tun¬ 
nelling through a single-level quantum dot. This system is 
also known as a single-electron transistor at the nanoscale 
and has been discussed in Section IV B. The experiment in 
Ref. [26] has been carried out in the kHZ frequency range, 
where a time-resolved measurement of the current trace is 
possible. Although all the statistics—including correlation 
functions of arbitrary order as well as the WTD—can directly 
be computed from this time-resolved current trace, this ex¬ 
periment provides an ideal test-bed for our proposal. We can 
compare the WTD obtained from our reconstruction scheme 
based on cMPS with the WTD directly deduced from the ex¬ 
perimental current trace. 

Due to the simplicity of the setup, our proposed method 
to access the WTD only requires the two-point function C 2 . 
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t[s] 


FIG. 3. (Color online). Two-point correlation function C 2 for sin¬ 
gle fermions tunnelling through a single-level quantum dot with data 
from Ref. [26] (red-dotted curve). The blue curve, obtained from the 
reconstructed values of the parameters rL,R using a cMPS formula¬ 
tion of the quantum experiment, agrees well with the experimental 
measurement of C 2 . The deviation for small times is due to experi¬ 
mental limitations in the time bin with respect to which the current 
trace is resolved. 


This one can directly be derived from the experimental spike 
train / (the time-resolved current trace) and is shown in Fig. 3 
(red dots). The rates Fl = 13.23 kHz and Fr = 4.81 kHz 
have been determined experimentally and the corresponding 
C 2 -function agrees very well with the analytical expression 
when the detector rate is taken into account[26] 

In our reconstruction scheme, the quantity Fl + Fr can be 
determined from the current spike train autocorrelation func¬ 
tion / ^ / by least squares methods or spectral estimation pro¬ 
cedures analogous to the procedure described in Ref. [10]. By 
requiring Fl > Fr and using the expression of the steady- 
state current (see Eq. (30)), Fl and Fr can be uniquely iden¬ 
tified. The reconstructed values for the rates are 

F^con ^ 10.80kHz, (35) 

Fgcon_ 4 76kHz. (36) 

The differences to the values from Ref. [26] are well within 
the range we would expect, regarding the time-resolution in 
the spike train data. The curve plotted from these recon¬ 
structed values of the parameters Fl,r is shown in Fig. 3 in 
blue. The slight deviation between the experimental points 
and this reconstructed C 2 -function is due to the discretisation 
of the counting time intervals used in the experiment: The 
size of each time bin is not much smaller than the time scale 
on which C 2 changes mostly. This leads to an error in the es¬ 
timation of the damping factor Fl + Fr and explains the dif¬ 
ference of the blue and the red dotted curves. Naturally, one 
could expect a more accurate reconstruction of the parameters 
FL and Fr when increasing the time resolution of the current 
trace or of the measurement of C 2 . 

From Fl and Fr, the corresponding cMPS matrices i?L and 
i?R can be constructed, as well as the matrices M and D. In 



r[s] 


FIG. 4. (Color online). WTD obtained from state-of-the-art exper¬ 
imental measurements with data from Ref. [26]. The reconstructed 
WTD using Eq. (33) is shown in blue. It matches the WTD obtained 
directly from the time resolved experimental current trace (red dots) 
well. The deviation is due to the finite-sized time bin corresponding 
to the resolution of the current trace. A more accurate reconstruc¬ 
tion of the WTD is expected by increasing the time resolution of the 
current trace or of the measurement of C 2 . 


this simple case, we did not need to employ the whole recon¬ 
struction procedure from Ref. [10]. Indeed, it is clear from 
Eq. (34) that only two out of the four parameters that charac¬ 
terise the system appear: Cn only depends on the tunnelling 
rates Fl and Fr —the eigenenergies 0 and 5 of H^ys do not 
contribute [43]. This will in general not be the case. 

The matrices i?L and i?R give access to the matrices D and 
V by direct computation. Inserting the latter into Eq. (33), the 
WTD can be reconstructed and the result is plotted in Fig. 4 
(blue curve). In order to build confidence in our procedure, we 
compare this result with the experimentally-accessible WTD 
(red dots). Let us recall that the transport rate is in the 
kHz range, hence the WTD can directly be extracted from 
the current spike train /: By sorting, counting all (discrete) 
waiting times between two consecutive incidents, and subse¬ 
quently normalising the resulting histogram, one obtains the 
red-dotted WTD in Fig. 4. The slight deviation between the 
WTD reconstructed via our proposal and the experimental one 
is again due to the discretisation of the counting time intervals. 
One could expect a more accurate reconstruction of the WTD 
when increasing the time resolution of the current trace. 

Although the WTD shown in Fig. 4 refiects the most el¬ 
ementary transport properties of single independent particles 
through the Poissonian distribution, it bridges the gap between 
theoretical predictions and experiments. It demonstrates that 
our reconstruction procedure based on a cMPS formulation of 
an open quantum fermionic system is reliable to access the 
WTD from the measurements of low-order correlation func¬ 
tions. This opens the route to access the WTDs in the quantum 
regime from low-order correlation-functions measurements in 
the high-frequency domain. 
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VI. CONCLUSION 

In this work, we have taken an approach motivated by 
continuous matrix product states to perform tomographic re¬ 
constructions of quantum transport experiments. On a for¬ 
mal level, we have extended this formalism to perform a 
reconstruction of unknown dissipative processes based on 
the knowledge of low-order counting probabilities. We then 
demonstrated that continuous matrix product states is an ad¬ 
equate formalism to describe quantum transport experiments 
based on tight-binding Hamiltonians. 

This work advocates a paradigm change in the analysis of 
transport experiments. The traditional method is to make ex¬ 
plicit use of a model to put the estimated quantities into con¬ 
text, a model that may or may not precisely reflect the physi¬ 
cal situation at hand. The cMPS approach is to not assume the 
form of the model, with the exception that the quantum state 
can be described by a cMPS. Such an approach is of particu¬ 
lar interest as it opens the way to the access of quantities that 
are not measurable experimentally with current technologies, 
high-order correlation functions and distributions of waiting 
times. 

To convincingly demonstrate the functioning of cMPS to¬ 
mographic tools applied to quantum transport experiments, 
we presented a simple example that consists of electrons tun¬ 
nelling through a single-level quantum dot. Making use of ex¬ 
perimental data, we showed that we could successfully recon¬ 
struct the distribution of waiting times from the measurement 
of the two-point correlation function only. This work consti¬ 
tutes therefore a significant step towards accessing the wait¬ 
ing time distribution in the quantum regime experimentally, a 
challenge present for several years now. Importantly, the ap¬ 
plication of our reconstruction procedure goes beyond the in¬ 
terest in waiting time distributions: It also provides an access 
to higher-order correlation functions, which are key quantities 
to better understand interacting quantum systems. 

In subsequent research, it would be desirable to further flesh 
out the statistical aspects of the problem. After all, the de¬ 
scription in terms of continuous matrix product states con¬ 
stitutes a statistical model. It would constitute an exciting 
enterprise in its own right to identify region estimators that 
provide efficiently computable and reliable confidence regions 
[44] when considering the problem as a statistical estimation 
problem, related to the framework put forth in Refs. [45-47]. 
We hope that the present work inspires such further studies of 
transport problems in the mindset of quantum tomography. 
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Appendix A: Reconstruction method from low-order counting 
probabilities 

In this appendix, we provide further technical details on 
the reconstruction scheme based on the measurement of low- 
order counting probabilities, Pq and Pi . The goal is to access 
the central cMPS parameter matrices M. and V. We refer to 
Fig. 1 for a general view of the reconstructible items. 

We start from Eq. (13) in the main text. By changing the 
integration bounds, we obtain Eq. (14), 

Pnir) = / dil / dt2--- dtnCn{T,ti,t2,...,tn), 

Jo Jo Jo 

where the integrand Cn is altered to 

CniT,ti,...,tn) ■■= lim R„--- 

L^oo 

■■■Rl® (g) . (Al) 

Note that in contrast to Eq. (7) where the propagating matrix 
is the transfer matrix T defined by Eq. (2) (or equivalently 
its diagonal representation D), the propagating matrix in the 
exponential terms between two measurement points now is the 
matrix S, which is defined by 

S:=Q* ®1 + 1®Q = T -^R*®Rj. (A2) 

j 

We can further simplify Eq. (Al) by performing the thermo¬ 
dynamic limit L ^ oo. The spectrum of T for a generic 
system consists of complex values with negative real part and 
only one eigenvalue being equal to zero. When taking the 
limit L ^ oc, all eigenvalue contributions to vanish, 

except the one corresponding to the zero eigenvalue. Hence 

lim = XeielX -^, (A3) 

L —>-CXD 

with the first canonical unit vector denoted by ei and the basis 
transformation matrix X to the diagonal basis of T. Similar 
to D and T, we define the matrix V as the diagonal matrix of 
S with basis transformation matrix Y 

S = YVY -^. (A4) 

By defining the matrices J\Aj := (g) Rj)Y for j = 

l,...,n, and setting Z := {zj^k) •= Y~^X (with inverse 
=: (%,/c)), we arrive at Eq. (15), 

Cn {pPl'i ' ' ' Pn) — 

Mie^C-ti)zei. 

Eq. (15) and Eq. (7) have a close structural resemblance: the 
matrices M and A4 are similar in the linear algebra sense, 
i.e., there exists a basis transformation from JVI to M. The 
matrices D and V are the diagonal matrices of the transfer 
matrix T and the matrix S, respectively. It is straightforward 
to transform M and D into JYl and V and vice versa: By 
subtracting M from D, we obtain S (up to similarity/basis 
transformation), whose diagonal matrix is V. Applying the 
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same basis transformation (from [D — M)ioV)io the matrix 
M results in the matrix Af. 

For n = 0, the counting probability function then reads 

Po(t) = elX-^Ye^^^Y-^Xei , (A5) 

which can be rewritten as a sum of complex exponential terms, 
with {fij} being the eigenvalues of S as 

Po{t) = • (A6) 

i=i 

This expression corresponds to the analogue of Eq. (10) in the 
main text. Since S is by definition a Kronecker sum of Q* and 
Q with eigenvalues {g*} and {g^} respectively, the spectrum 
of S consists of the sums g* + qk with j^k = 1,..., d. It is 
closed under complex conjugation (for each element of the set 
its complex conjugate is also element of the set), as well as the 
coefficient set {zijZj^i}. This ensures that Pq is real-valued. 
Being related to Q (which consists of a skew-hermitian matrix 
(with imaginary spectrum) and negative definite matrices), we 
have that Re jij < 0 for each j, such that all summands van¬ 
ish sufficiently fast and Pq is normalisable. Furthermore, the 
dominance of the damping factors over the oscillatory compo¬ 
nents ensures the positivity of Pq (in particular, the fij with the 
least damping is always real-valued). Analogously, for Pi(r) 
we obtain 


j,k^k,l 

j,k=l 


(1 ^j,k) 


0/^fcP _ 

k'k 


with the Kronecker delta 





(A7) 


argument: All probability functions Pn are invariant under 
scaling and permutation of the eigenvectors in the matrices 
X and Y (except for the eigenvector of T corresponding to 
eigenvalue zero). This allows us to require all but one zi^j to 
be equal to one, and immediately obtain the according number 
Zj^i. The remaining coefficient can then be determined via the 
normalisation constraint 


P 

(A9) 

i=i 

so that all Zj^i are known. This can be used to obtain the 
diagonal elements Mj^j from Af. For the remaining matrix 
elements, only the symmetric elements M.j^kPM.k,j (but not 
their constituents) are directly accessible since 


jAf j, 


kZ'k,l~ 


j,k=l 


k'k 


/ (yZijXij^k^k,! P /cj’2^jf^i) . (AlO) 


j<k 


k'k k'j 


However, this does not constitute a limitation for the recon¬ 
struction of the matrices Q and R of the ancillary system. To 
this end, we make use of the inner structure of Af. The diago¬ 
nal matrix V with eigenvalues fij can be reordered such that it 
has the form V = Dq (g) 1 -]- 1 (g) Dq with diagonal Dq con¬ 
sisting of the eigenvalues of Q. Reordering the eigenvectors 
in Y accordingly, we can assume that the matrix Y and hence 
the matrix A1 have the form of a Kronecker product 


A1 = ^ ^rec . (All) 


/('t) = 

m,n 

Assuming that the terms are linearly independent, 

in principle one can always single out these contributions as 
well as their corresponding prefactors Cm,n- This gives us the 
chance to extract the coefficients {zijZj^i} and the eigenval¬ 
ues {fii} from Pq, provided that no coefficient is identical to 
zero. Rearranging the values {/i^} to a diagonal matrix in Kro¬ 
necker sum form results in the matrix V. One should note, 
however, that efficient spectral recovery algorithms like the 
matrix pencil method do not straightforwardly work for func¬ 
tions such as Pn, n > 2, where the exponential functions are 
multiplied with powers of r. 

In order to reconstruct the elements of the matrix A1 to¬ 
gether with the off-diagonal elements of Z, we use a gauge 


Here Prec = ^ ([^dxd general not diagonal. 

The symmetrised components of Af can then be written as 
and the constituents rj can be determined 
(up to a phase factor) by equating them with the coefficients in 
Eq. (AlO). The according equation system can then be solved. 

The important point is that Prec and Qrec •= Dq are valid 
cMPS parameter matrices in the same gauge and hence are 
sufficient for reconstruction with the same argument as in Ref. 
[10, III.E]. Let us note that concrete values of the basis trans¬ 
formation matrices X and Y are in fact never used or needed 
in the reconstruction procedure. From Prec and Qrec. we can 
compute all quantities we need to establish the correlation and 
counting probability functions, in particular Af and V. Re¬ 
gauging Prec and Qrec such that the orthonormalisation con¬ 
dition [1] is fulfilled, yields a reconstruction of the free Hamil¬ 
tonian Prec of the ancillary system. 
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